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We argue that both the positron fraction measured by PAMELA and the peculiar spectral features reported 
in the total electron-positron flux measured by ATIC have a very natural explanation in electron-positron 



■ pairs produced by nearby pulsars. While this possibility was pointed out a long time ago, the greatly 

q , improved quality of current data potentially allow to reverse-engineer the problem: given the regions of 

in : 

c/5 pulsar parameter space favored by PAMELA and by ATIC, are there known pulsars that explain the data 

with reasonable assumptions on the injected electron-positron pairs? In the context of simple benchmark 
models for estimating the electron-positron output, we consider all known pulsars, as listed in the most 

t> . 

complete available catalogue. We find that it is unlikely that a single pulsar be responsible for both the 
PAMELA positron fraction anomaly and for the ATIC excess, although two single sources are in principle 
enough to explain both experimental results. The PAMELA excess positrons likely come from a set of 

oo ; 

, mature pulsars (age ~ xlO yr), with a distance of 0.8-1 kpc, or from a single, younger and closer source 

like Gcminga. The ATIC data require a larger (and less plausible) energy output, and favor an origin 

"X' 

associated to powerful, more distant (1-2 kpc) and younger (age ~ 5 x 10 yr) pulsars. We list several 
candidate pulsars that can individually or coherently contribute to explain the PAMELA and ATIC data. 
Although generally suppressed, we find that the contribution of pulsars more distant than 1-2 kpc could 
contribute for the ATIC excess. Finally, we stress the multi-faceted and decisive role that Fermi-LAT will 
play in the very near future by (1) providing us with an exquisite measurement of the electron-positron flux, 
(2) unveiling the existence of as yet undetected gamma-ray pulsars, and (3) searching for anisotropies in the 
arrival direction of high-energy electrons and positrons. 

PACS numbers: 97.60. Gb, 98.70.Sa,95.85.Ry 
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I. INTRODUCTION 



'positrons over the sum of 
l| and of the electron-plus- 
3[ stirred 



Recent measurements of the positron fraction (the fraction o 
positrons and electrons) reported by the PAMELA collaboration 
positron differential energy spectrum reported by ATIC and by PPB-BETS 
up great interest both in the cosmic-ray community and among those who work in the field 
of indirect dark matter searches. 

The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics 
(PAMELA) measured a statistically significant trend of increasing positron fraction with 
energy in the 10 < E e ±/GeV < 100 GeV range, at odds with the standard expectation in 
the context of secondary positron production from interactions of cosmic-ray nuclei with the 
interstellar gas: if most positrons are of secondary origin, the positron fraction should fall 
as a function of increasing energy (see e.g. the recent comprehensive stuy of Ref. Q]). In the 
past, a similar trend had actually been repeatedly observed, although with much lower sta- 
tistical significance, by other experiments, starting as early as 1969 More recent results 
pointing towards an increasing positron fraction above 10 GeV include the measurement 
reported in Ref. 6] (1975), Ref. Q, Ref. 8] (1994), the data reported by HEAT (1997, 
Ref. 9|n and those quoted by the Alpha Magnetic Spectrometer (AMS) collaboration (2000, 
Ref. (lo|). 

The Advanced Thin Ionization Calorimeter (ATIC) recently reported a "bump" in the 
high energy flux of electrons and positrons in excess of the standard expectation, and 
marginally in agreement with a statistically less significant excess also visible in the PPB- 
BETS data reported in Ref. [sj] 1 . Outside the bump, the ATIC data are in substantial 
agreement with previous measurements, including those reported 
older PPB-BETS data [l^, as well as emulsion chambers data 
spectral feature observed in the data from both ATIC flights, and which had never been 
reported before, is a sharp excess in the electron-positron flux, specifically in the energy 
range between 300 and 800 GeV 

In addition to the mentioned experimental results on cosmic-ray electrons and positrons, 
the atmospheric Cherenkov telescope (ACT) H.E.S.S. recently reported measurements of 



by AMS [lfj], BETS [llj, 
la . \l4\. The anomalous 



PPB-BETS measured 84 electron events above 100 GeV Of, while ATIC 1724 



a, 
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the flux of very high-energy (4£e± > 600 GeV) electrons, with data extending all the way 
up to energies around 5 TeV [15|. The data were taken during off-target observations of 
extragalactic objects. Since no conclusive distinction is possible with ACT's between e 
and gamma-rays, the data likely include a contamination from diffuse extra-galactic (and 
possibly, to some extent, diffuse galactic) gamma rays as well as from residual hadronic 
cosmic-ray events. The H.E.S.S. data should therefore properly be regarded as upper limits 
to the e ± flux at large energies, and this is how they will be considered in the present study. 

The nature of cosmic-ray positrons have long been believed to be dominantly of sec- 
ondary origin: the observed cosmic-ray positrons, in that scenario, are the yield of inelastic 
interactions of (primary and secondary) cosmic-ray protons with the interstellar gas (see e.g. 
Ref. 16], 17] and [3]). Calculations of the cosmic-ray positron and electron flux as well as 



of the fraction of secondary positrons, assuming several propagation models and diffusion 
setups, have been carried out extensively over the last 30 years (see e.g. [lj]], [^(J, [l]), with 
increasing degrees of sophistication, and relying on experimental data of better and better 
quality . The generic prediction of all these models is a slowly decreasing positron fraction 
with increasing energy, and of a steeply declining high energy differential primary electron 
flux (for instance, the conventional and optimized models of Strong et al. 21|, which give a 
diffuse galactic continuum gamma-ray flux compatible with EGRET data as well as cosmic- 
ray abundances compatible with a wide set of data, feature injection spectral indexes for 
high-energy electrons of 2.54 and 2.42, respectively). 

As mentioned above, however, experimental data on cosmic-ray positrons have actually 
shown for quite a long time a discrepancy with a purely secondary origin for positron cosmic 
rays. For instance, Ref. 22J pointed out almost 20 years ago that, in addition to diffuse 
sources such as secondary cosmic-ray production and to primary sources, such as electrons 
re-accelerated in supernova remnant (SNR) shocks, data on the positron fraction and on 
the flux rather conclusively pointed towards the existence of a powerful, nearby source 
of energetic primary electrons and positrons above 20 GeV. 



As early as 1995, and commenting on data from before 1990, Ref. [23J stated that "the 
measured content of positrons in the total electron flux, [. . .], at least at high energies [. . .], 
is a[n] [...] enigma awaiting an explanation. [...] it is obvious that some other source of 



positrons is needed 1 . Ref. 



23j also argued for the first time that nearby pulsars might seed 



a high-energy positron excess, as well as non-standard spectral features in the cosmic-ray 
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electron-positron spectrum. The positron anomaly is therefore nothing qualitatively new. 
We argue here, however, that the greatly improved quality of recent data allows us to 
quantitatively push previous analyses to the level of associating electron-positron data to 
existing, observationally well-known galactic objects. This is all the more exciting given the 
apparently anomalous data reported in the total differential electron-positron flux. 

A physical hallmark of energetic (E e ± > 10 GeV) electrons and positrons is that they 
loose energy very efficiently via two dominant processes: inverse Compton (IC) scattering 
of background radiation photons, and synchrotron radiative losses. The energy loss rate for 
both processes goes as the square of the e energy. Defining the lifetime t of relativistic 
e* 1 in the interstellar medium as the time it would take for them to cool from an injection 
energy E e ± to rest energy (E e ± = m e ) if they only lost energy via IC and synchrotron losses, 
one finds 

t ~ 5 x 10 5 ( — ^ yr. (1) 



E e ± J 

This implies that TeV e^ 1 must have been injected not much longer than ~ 10 5 yr ago. This 
also implies that, for conventional values of the diffusion coefficient D for the propagation of 
electrons and positrons, they must have been produced somewhere within \/D x t ~100-500 
pc. The efficient energy losses of energetic electrons therefore constrain the sources of high 



energy cosmic-ray electrons and positrons both in time and in space [24J . Any interpretation 
of the reported "lepton anomalies" must evidently be compatible with this fundamental 
physical property of cosmic-ray electron and positron propagation in the interstellar medium 
(ISM). 

Several scenarios have been suggested to explain a deviation of the electron-positron data 
from the standard cosmic ray secondary positron model (for a recent model-independent 



analysis see e.g. Ref. 



25]). A list of these includes: 



1. A cutoff in the primary electron spectrum at large energies (see e.g. 26fl). This 
possibility might help explain the positron fraction data, but is evidently not an ex- 
planation to the measured differential flux, which would still be dominated, at high 
energy, by primary cosmic-ray electrons; 

2. The production and re- acceleration of secondary electrons and positrons via hadronic 
processes in giant molecular clouds. For instance, for nominal model parameters, 



Ref. 



271 ] points out that those systems can source a significantly larger positron frac- 



tion in excess to the prediction of ordinary diffuse cosmic-ray models, starting at 



energies around 10 GeV. In this scenario, the results of Ref. [28| indicate that it could 
be possible to reproduce the PAMELA data below 30 GeV, but no satisfactory expla- 
nation to ATIC or to the higher energy bins reported in the PAMELA data seems to 
be achievable in the context of giant molecular clouds models; 



3. Ref. 



29( pointed out that might be pair-produced near compact gamma-ray sources, 



for instance as a result of processes involving, in the initial state, the source gamma 
rays and the optical and UV photons also produced locally by the gamma-ray emitter. 
Electron-positron pair production could also occur from compact gamma-ray sources 
in interactions of gamma rays with starlight photons in the interstellar medium, as 
observed in Ref. 30[. This possibility appears, however, to be tightly constrained by 
current galactic gamma-ray data, and it will be further constrained and quite easily 
scrutinized by upcoming Fermi data. 

4. Another possible source of primary positrons is the radioactive f3 + decay of nuclei (such 
as 56 Co) ejected during a supernova explosion 31_|, |32j]. Supernova models, however, 
indicate that the spectrum of such nuclei (and therefore of the positrons produced 
in the decays) is cut off at energies around 100 GeV, ruling out this scenario as an 
explanation to the spectrum measured by ATIC. Also, the energetics of the process 
indicate that this mechanism is likely insufficient to produce as large a flux of positrons 
as to explain the PAMELA positron fraction data too. 

5. Pulsars are undisputed sources of electron-positron pairs, produced in the neutron star 
magnetosphere and, possibly, re-accelerated by the pulsar wind and/or in SNR shocks. 
Nearby pulsars have since long been considered as very likely dominant contributors 
to the local flux of high energy electrons and positrons (see e.g. Ref. 33|, |34|). As 
shown in Ref. 35J and 36], the recent PAMELA data can rather well be interpreted 



as originating from one single nearby pulsar, e. g. G eminga, or by the co 
position of the e ± flux from all galactic pulsars [3^, (see also Ref. 



lerent super- 
for a recent 



re-assessment of pulsars as the origin of the cosmic-ray lepton anomalies appeared 
after the present manuscript was first posted); 



6. A further possibility that recently received enormous attention is that the excess 
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positrons and electrons originate from the annihilation or decay of particle dark mat- 



ter (for a possibly incomplete list of related studies appeared 



aefore the first ver 



sion o 
47 



68 



' the present manuscript was released see Ref. |39 
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the present v ersion of this manuscript include 
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123, 



137|]). As opposed to 



all of the other possibilities mentioned above, a dark matter interpretation invokes an 
entity whose fundamental particle physics nature has yet to be unveiled, and whose 
existence in the form of a particle with a mass in the GeV-TeV range is yet to be 
demonstrated. In addition to this, there are at least three general weaknesses to the 
dark matter interpretation of the positron fraction and of the features detected in the 
high energy electron-positron flux: 



(a) the annihilation or decay products of dark matter generically include both leptons 
and hadrons. Unless the dark matter model is tuned in such a way as to make it 

n 

"lepto-philic" [66|, the hadronic products will typically overpro duce both gamma 
rays and antiprotons, in contrast with experimental data (e.g. 138| ); 

(b) the annihilation rate required to explain the cosmic-ray data under discussion is 
typically a couple orders of magnitude larger than what expected if the dark mat- 
ter is cosmo logically produced through thermal freeze-out in the right amount. 
This implies a further model-building intricacy. In addition, low-velocity en- 
hancements, which would lead to a burst of annihilations in the first collapsed 

1391 ]. are generically ve ry tig htly constrained by observations of the cosmic 



halos 



diffuse background radiation 140 ] 



(c) The large amount of energetic positrons and electrons injected in our own Galaxy 
and in any dark matter structure in the Universe would yield secon dary emis- 
sion at all wave-lengths (for a recent comprehensive analysis see Ref. 1411 ]): this 
includes radio signals from synchrotron losses, X-rays and gamma-rays from in- 
verse Compton scattering off the cosmic microwave background and other back- 
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ground radiation photons, as well as fr om b remsstrahlung, neutral pion decay 



and internal bremsstrahlung radiation |142j |. Tight constraints exist at each 
one of these wavelengths, at all dark matter halo scales and distances, rang- 
ing from local dwarf galaxies to distant clusters (a very incomplete list of refer- 



ences on mult i- wave 
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144 
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145 
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146. 
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i matter annihilation inc 
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Recent studies, appeared after this manuscript was first released, discuss further possib 
explanations to the cosmic ray "lepton anomalies" . These include: a gamma-ray burst 164 ] 



inelastic collisions between high energy cosmic rays and bac kground radiation photons in 



the astrophysical sources where cosmic rays are acceler ated |165l . 11661 ] . inhomogeneities in 
the spatial distribution of sources of cosmic-ray leptons 1671 ] , hadronic processes prod ucing 
seco ndary positrons in the very same sites where primary cosmic rays are accelerated 168 . 
1691 ] and s scenario where one (or more) recent suernova event (s) occurred in gas clouds 
close to the Solar System [170 ]. 

Of all possibilities listed above, we focus in the present study on the one which we believe - 
and intend to argue with the present analysis - is the best motivated: nearby, known galactic 
pulsars. On the one hand, pulsars are predicted to produce e ± pairs with a spectrum that 
is suitable in principle to explain both the data on the positron fraction and, possibly, the 



overall electron-positron differential flux. On the other han d, pu 



well studied objects, and catalogues of galactic pulsars exist 



171 



sars are well known and 



1721 1 , off of which we can 
compare theoretical speculations on the most plausible sources for the observed e* 1 fluxes 
with actual observational data. 

Unlike in the dark matter hypothesis, the pulsar scenario does not require the introduction 
of an extra, exotic component to explain the data. This goes well with the celebrated quote 
of the 14th-century English logician and Franciscan friar, William of Ockham: ll Entia non 
sunt multiplicanda praeter necessitatem" : there is no need to introduce extra entities beyond 
necessity. 

The scope of the present analysis is to "dissect" the PAMELA and ATIC data in the spirit 
of (or with) Occam's razor: we will pinpoint which regions of the pulsar parameter space 
(to be defined in the next section) are favored by experimental data; we will then turn to 
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observations, and see if those parameter space regions correspond to existing astrophysical 
objects, and if those objects are predicted to have an e ± energy output in line with what is 
needed to explain the cosmic ray data. We shall argue that the recent positron and electron 
data are certainly very intriguing, but that they are not anomalous, since a very natural 
interpretation exists in terms of known, existing astrophysical objects. We will also outline 
the major uncertainties that enter the process of reverse-engineering cosmic-ray lepton data 
to infer their possible origin. Among these, uncertainties in the pulsars' age and location 
at the point in time when the bulk of the high-energy lepton pairs where injected play an 
extremely relevant role. 

In summary, highlights and novelties presented in this study include: 



accounting for all known, existing pulsars, as listed in the ATNF catalogue [171] : 
this approach allows us to (i) understand, qualitatively, which regions of the pul- 
sar age-versus-distance parameter space contribute in given ranges of the cosmic ray 
electron-positron spectrum (see sec lVIIII and fig. fT2l and fT3l) . (ii) analyze, through the 
determination of the relevant age-distance range, where the pulsar catalogue incom- 
pleteness can affect predictions for the high-energy cosmic-ray electron-positron spec- 
trum, and thus assess where and how newly discovered gamma-ray pulsars can impact 
our understanding of the origin of high-energy cosmic-ray electrons (see e.g. sec JIXI and 
fig. [15]) . (iii) address how the hereby proposed "reverse-engineering" approach (and 
in general any attempt at pinpointing a local high-energy electron-positron source as 
the main contributor to the leptonic cosmic-ray fluxes measured at Earth) is affected 
by uncertainties in the determination of the pulsars' age and distance, due to e.g. to 
trapping of the in the pulsar wind nebula, and to pulsar kick velocities (see fig. [TOT) . 
(iv) estimate, by considering all ATNF pulsars, the theoretical uncertainties built in 
several current "benchmark" models for the e ± energy injection from pulsars, such as 
those described in the next section (see fig. [Hand EJ); 

systematically analyzing and quantifying the role of uncertainties related to cosmic- 
ray energy losses and diffusion modeling (see e.g. figlMH]), and quantifying how these 
uncertainties impact the identification of pulsars that can source the positron fraction 
data and/or the spectral features; 

pointing out that there are regions of the pulsar age-distance parameter space that are 



more likely to contribute to the detected lepton anomalies (see e.g. figj9]and tab. IIIII) . 
This point is relevant to the possibility of carrying out anisotropy studies, since it 
provides e.g. target directions for the possible detection of a dipolar asymmetry in 
the CRE arrival directions. In addition, as we point out with figHH Fermi data may 
soon tell us whether one pulsar or a combination of pulsars are at the origin of the 
high-energy end of the cosmic-ray electron-positron spectrum, and has the potential 
to actually detect the mentioned anisotropics; 

• theoretically studying the role of leaky-box boundary conditions to the widely-used 
analytical expression 23] employed to calculate the spectrum, at Earth, of cosmic- 
ray electrons and positrons injected by a pulsar in the burst-like approximation; this 
physically means that we provide an estimate of the fraction of electrons and positrons 
injected by a given pulsar at a given point in time that freely escapes the galactic 
diffusive region, as a function of the pulsar and of the diffusion model parameters (see 
the last few paragraphs of sec. [Til and fig. [3]). 

It goes without saying that the quest for the fundamental nature of particle dark matter, 
in this case indirect detection, needs to undergo the perhaps painful, or prosaic, process of 
deep scrutiny of any "mundane", standard interpretation of experimental data, including 
known astrophysical processes. The present study is a modest attempt in this direction. An 
important point we wish make at the end of this study is that the Fermi gamma-ray space 



telescope will play a decisive role, in the very near future, to elucidate the natu re o 



the e d 



sources, as well as to inform us on what particle dark matter can or cannot be [1731 ] . 

The outline of this study is as follows: in the next section we describe the processes 
that can seed pair production in pulsars, estimate the e ± energy output and injection 
spectrum, and discuss the subsequent particle propagation, including a novel analytical 
formula to account for leaky-box boundary conditions in the propagation setup. Sec. IIIII 
gives a few examples of nearby pulsars that could contribute to the flux observed locally 
at the level of the measurements reported by ATIC and PAMELA. The following section 
addresses the question of whether those measurements can be explained by a single nearby 
source. Sec. [V] and [VTJ are then devoted to outlining the pulsar parameter space favored 
by the PAMELA and ATIC data, respectively. Sec. IVHI contrasts a scenario where the 
measured e* 1 originate from a single source to one where they stem from multiple pulsars in 
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well defined parameter space regions. We carry out a self-consistency check of our models 
in sec. IVIIIt where we compute the e summing over the full galactic pulsar catalogue. 
Finally, we outline in sec. IIXI the decisive role of the Fermi telescope in understanding the 
data: (1) to provide an exquisitely detailed measurement of the e^ 1 spectrum from 20 
GeV to 1 TeV, and (2) to discover potentially relevant new gamma-ray pulsars, and , and 
(3) to search for anisotropics in the arrival direction of high-energy electrons and positrons. 
Sec. |X] summarizes and concludes. 



II. BENCHMARK MODELS FOR THE PRODUCTION OF ELECTRON- 
POSITRON PAIRS FROM PULSARS 



Pulsars are known to be potentially very powerful sources of primary electrons and 
positrons, pair-produced in the neutron star magnetosphere. Two classes of mechanisms of 
pulsar relativistic particle acceleration (and associated gamma-ray and e 1 * 1 pair-production) 
have been widely discussed in the literature, corresponding to the regions of the mag- 



netos pher e where the p a rtic 



Ref. 



174 



175 



176 



177 



178 



e acceleration takes place: the volar cay mod el (see e.g 



179| ) and the outer gap model 



180 



181 



1821). We will 



not review those mechanisms in any detail here. In short, in both models very high energy 
electrons are accelerated by electric fields in the pulsar magnetosphere. The electrons, in 
turn, synchrotron-radiate gamma rays with energies large enough to pair-produce electron- 
positron pairs in the intense magnetic fields of the magnetosphere, and/or through scattering 
off of intervening photon fields, such as the thermal X-rays produced by the pulsar itself. 
We only mention here that, in the context of the outer gap model, the threshold for e ± 
pair-production from the scattering off of the pulsar thermal X-rays onsets when the ratio 
of the siz e of the ou ter g ap versus the radius of the light cylinder is less than one (for details 



see Ref. 



1831 ] and 



184J). Pulsars that fulfill this condition are indicated as "gamma-ray 



pulsars". In terms of the pulsar magnetic field at the neu tron star surface £>i2, in units 
of 10 12 G, and of the pulsar period P in seconds, Ref. 
correspond to those objects that fulfill the condition 



184j finds that gamma-ray pulsars 



g = 5.5 P 26/21 B- 2 i/7 < 1. 



(2) 



A crucial parameter for the present study is the estimate of the total energy that a pulsar 
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radiates in electron-positron pairs that actually diffuse in the interstellar medium (ISM). To 
estimate the total energy output in e from pulsars with rotation period P/s (corresponding 
to a rotational frequency fl = 2ir/P Hz) and a period d eriva tive P , we adopt her e th e 



classic magnetic dipole pulsar model outlined e.g. in Ref. [185| , and reviewed e.g. in 186 ] . 
In this scheme, pulsars consist of a rotating neutron star featuring a dipolar magnetic field 
which is not aligned with the rotation axis. The neutron star, with a moment of inertia 
/ ~ 1.0 x 10 45 g cm 2 and angular frequency fl is assumed to lose rotational energy via 
magnetic dipole radiation, and to follow a braking law of fl oc — fl 3 . This can be integrated 
to yield (in the limit of a mature pulsar, i.e. one with a rotation period much longer than 
its initial one) a characteristic age T — ft/ (2ft). 

Notice that this estimate of the pulsar age does not necessarily correspond to the actual 
pulsar age (some of the pulsar ages in the ATNF catalogue obtained in this way, in fact, turn 
out to be even larger than the age of the Universe!). T his p roble m is particularly severe for 
young pulsars with ages between 10 4 and 10 5 yr. Ref. [187| and [1881 ] . for instance, pointed 
out that the pulsar J181 1-1925 associated to the SNR Gil. 2-0. 3, whose age is inferred from 
the likely association with the supernova reported by Chinese astronomers in A.D. 386, has 
an actual age which is more than one order of magnitude smaller than the characteristic 
age derived from timing measurements, w hich is of the order of 24,000 yr. From a semi- 



analytical standpoint, in the model of Ref. 



1891 ] the median over-estimate of the pulsar age 



obtained from timing measurements is of 60,000 yr, which gives a yard-stick to estimate 
whether for a given set of pulsars this uncertainty is or not relevant. 

In addition to the uncertainty in the actual pulsar age, another effect should be taken 
into account when considering a burst-like injection of e ± from pulsars: before diffusing 
in the ISM, the electrons and positrons produced by the pulsar are trapped in the pulsar 
wind nebula or in the supernova remnant that envelopes it. This time-lag is of th e order 
of 10 4 ...10 5 yr (for a review on pulsar wind nebulae we refer the reader to Ref. 190I ]). This 
time-delay also implies that pulsars relevant for the present discussion are mature pulsars, 
with typical ages in excess of 10 5 yr - younger pulsars likely having not expelled yet the 
produced electron-positron pairs in the ISM. 

We outline below here four possible schemes to estimate the total energy injected in 
the ISM by pulsars in the form of electron-positron pairs, and evaluate the corresponding 
"theoretical uncertainty" by evaluating the output of each model for all pulsars listed in the 
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most complete available catalogue. 

The spin-down power of a pulsar corresponds to the pulsar's rotational energy dissipated 
in dipolar electro-magentic radiation, and is given by E = IQti. The dipole m agne tic field 



strength at the stellar surface, again in the classic pulsar dipole model of Ref. 



1851 ] . reads 



B = 3.2 x 10 19 \AP x P G. (3) 

We indicate with f e ± the fraction of the rotational energy (integrated spin-down power) 
which is eventually dissipated in e ± pairs diffusing in the ISM. We can naively estimate 
the energy output as follows: The general solution to the pulsar rotational velocity, in the 
regime where the pulsar looses energy dominantly via dipole radiation, can be cast as 

Qo_ 
t/roY 

where Q is the initial spin frequency of the pulsar, and, indicating with R s the neutron star 
radius, 

3c 3 / 

is a characteristic "decay time" which, incidentally, cannot be directly derived from pulsar 
liming measurements, and is usually assumed to be, for nominal pulsar parameters, ~ 10 4 yr 



m = ~ .° a/ 2 . ( 4 ) 



24j . For known mature pulsars, which could contribute to the observe d flux because the 



produced particles are no longer trapped and lose energy in the nebula 1911 ] . we approximate 
the total energy output as 

£ out = / e ± x / x Jdtnn^ ^-m 2 (6) 

Now, since 

r 



fi^fi 2 -, (7) 



we have, for t — T, i.e. the age of the pulsar, 

E out [ST] = ^m 2 - = f e± (im (§:) - = fe±E—. (8) 

In this naive scenario, the pulsar output in e* 1 is then completely fixed by the rotation 
period and its derivative through the derived values of the spin down power E and typical 
age T, once the efficiency factor for e ± production is taken into account. We indicate 
hereafter this model as the "ST" model. Unless otherwise specified we assume f e ± = 0.03, 
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i.e. 3% of the rotational energy gets converted into ejected e ± pairs. This is a rather 



conservative assumption, along the lines of what assumed e.g. in Ref. 



351 ] . i.e. / e ± ~ few 



%. A quantitative discussion of plausible values for f e ± was recently given in Ref. 38] . 



We shall not review their discussion here, but Ref. [38| argues (see in particular their very 
informative App. B and C) that in the context of a standard model for the pulsar wind 
nebulae, a reasonable range for f e ± falls between 1% and 30%. 

An altern ative scheme is the model for the pulsar energy output formulated by Harding 
and Ramaty 1921 ] (hereafter, the HR model), where the main assumption is that the pulsar 
luminosity in e ± is proportional, via a factor fna, to the gamma-ray output (the high-energy 
gamma rays and e ± are assumed to follow the same spectral index): 



L e ± = Ihr x L 7 , where 
~ 7.1 x 10 33 (£/GeV)- Q B 12 P~ L7 ph s" 1 GeV^ 1 



(9) 



The expression given above is, incidentally, i n agr eement with data from the Crab and 
from the Vela pulsar, as pointed out in Ref. 122] . HR then assume that positrons are 
produced with gamma rays in the cascade processes described above. Particle Monte Carlo 
simulations indicate that in electro-magnetic showers the ratio of e pa irs t o gamma-rays 
fun lies between 0.2 and 0.5 in the energy range between 1 and 10 GeV 176| . Substituting 



the pulsar age (in s) and magnetic field for the rotation period (also in s) according to 

„ B 12 V2T 



3.2 x 10 



(10) 



and integrating over time, the total e energy output in the HR can be estimated as 

0.15 

erg. 



E out [HR] = 1.3 x 10 46 f HR Pf 2 l 



-0.7 



T 



10 4 yr 



(11) 



In the outer gap model of Chi, Cheng & Young (CCY) [1911 ]. the output can be 
assessed by taking the product of the total pair production rate in their setup (to which we 
refer the reader for additional details) 



3.4 x 10 37 g B™ /7 P" 8 / 21 s" 1 



(12) 



and of the typical electron injection energy 



l~l.lx 10 6 B% 7 P- 76 / 21 eV 



(13) 
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which, after substitution of Eq. (flOl) and after integration over time, gives: 

E out [CCY] = 7.6 x 10 48 C B; 2 9/7 (t^A ' erg, (14) 

where ( ~ 0.3 parameterizes the escape efficiency from the light cylinder. 

A similar but slightly different construction was proposed by Zhang and Cheng (ZC) 



371 ]: here, following their analysis and the procedure outlined above, we arrive at an energy 



output, for the ZC model, of 



E out [ZC] = 4.8 x 10 48 B~i (t^A erg. (15) 

We now proceed to assessing the e ± output for pulsa rs appearing in the Australia Tele- 
scope National Facility (ATNF) pulsar catalogue 2 171] ■ The ATNF sample provides the 
most exhaustive and updated list of known pulsars, including on only radio pulsars, but also 
those objects that were detected only at high energy (e.g. X-ray pulsars and soft gamma-ray 
repeaters). The total number of objects we use here (essentially all ATNF pulsars for which 
measurements of the pulsar period, its derivative and the distance are available) is 1789, of 
which 266 objects feature g < 1 ("gamma-ray pulsars" according to the definition above). 
The first point we wish to assess here is how large a range of predictions the four benchmark 
models outlined above give when applied to the entire ATNF pulsar catalogue. This will 
provide us with a rule of thumb on how uncertain predictions for the normalization of the 
contribution of known pulsars to the flux of local cosmic-ray electrons and positrons are. 

In fig. [1] we show the estimated pulsar energy output in for the entire ATNF sample, 
for the ST model (upper panels) and for the CCY model (lower panels). The two panels to 
the left indicate the pulsar position in the age versus e output plane (the energy output 
in ergs), while on the right we show the distance versus output plane. Red points highlight 
those pulsars whence a contribution from the outer gap is expected ( "gamma-ray" pulsars, 

g< !)■ 

The first feature we observe in our results is that for pulsars with ages 10 5 < T/yr < 10 6 
the model predictions fall between 10 47 and 10 48 erg. We remark here that these values are 



studies [35, 



36 



in agreement with older estimates, see e.g. 23|] and wit the estimates quoted in more recent 



38] . We notice that the relevant energy output range, quite remarkably, 



http : //www. atnf . csiro . au/research/pulsar/psrcat/ 



14 



Age [yr] 



H 

u 



10 



io J 



Distance [kpc] 

o 



3 10 



50 



■l-i 

3 

& 
3 

o 

M 

u 
e 
W 



10 
10 
10 
10 



49 



4X 



47 



T 1 — I I I I I I I 



ST model 



i — i — i i 1 1 1 1 



All ATNF pulsars 
g<l ("GR" Pulsars) 



• • • v "ft • 

o 



_i i i i ls 1 1 1 



_i i 



1 10 



50 



10 



49 



3 

o 

>, 

M 
u 
u 
5 
W 



10 
10 
10 
10 



47 



45 



CCY model 



_i i i 



10 




T 1 — I I I I I I I 1 1 — I I I I I I I 1 1 — I I I I lEj 



10 10 
Age [yr] 



10 

Distance [kpc] 



FIG. 1: Scatter plot in the energy output in e ± , in ergs, versus age (left panels) and versus distance 
(right panels), for the ST model and for the CCY model. The black circles indicate all pulsars in 
the ATNF catalogue, while the red dots indicate "gamma-ray" pulsars (g < 1). 

does not critically depend on the choice of the model for the estimate of the output. The 
spatial distribution of pulsars places the g < 1 outer gap candidates at around a kpc or so, 
with a larger scatter in the ST model (with young pulsars contributing lower outputs and 
older pulsars larger outputs). 

In the following fig. [2] we illustrate the ratio of the three HR, CCY and ZC model pre- 
dictions for the pulsar output over the ST prediction, as a function of age (left) and of 
distance (right). First, we notice that the HR model tends to predict a significantly lower 
output compared to other scenarios, as already noticed in the past • Secondly, compared 
to the CCY and ZC models, the ST model predicts a much larger output for mature pulsars 
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FIG. 2: The ratio of the predicted energy output from "gamma-ray" (g < 1) pulsars for the 
HR, ZC and CCY models over the naive prediction (ST model), as a function of age (left panel) 
and of distance (right panel). 

versus younger pulsars. In particular, it is easy to ascertain that the ratio of the ZC to the 
ST model scales exactly as T -1 , which quantifies the bias for the ST prediction towards 
older pulsars. No significant correlations with distance are observed for the ratios among 
the various models. 

Once the total output is specified, the injected electron-positron spectrum is fully de- 
fined by its spectral shape. It is well known that escaping the pulsar wind cannot be 
accelerated to arbitrarily high energies, and a cut-off is expected, depending on the pulsar 
environment, in the energy range around the TeV scale. The precise position of the cutoff is 
quite uncertain, and it also depend s critically on the pulsar age and on the magnetic fields 



in the pulsar wind nebula (see e.g. |l93l . Il94j ). Although physically there will be an upper 
limit to the energy of the escaping the pulsar, this would (i) introduce one more param- 
eter in the analysis we present below, and make it easier, in general, to accommodate the 
observed lepton anomalies by tuning it, and (ii) for the energies of interest and pulsars under 
consideration, the precise location of the cutoff is often not relevant (for instance because 
the pulsar is old enough that the maximal energy of propagating from it to the Earth is 
well below the injection spectrum cutoff). We therefore only illustrate the impact of a cutoff 
in the injection spectrum in a few selected cases below - if not specified otherwise, we will 



1(3 



not assume any injection spectrum cutoff in what follows. 

Below the spectral cutoff in the injected flux of e^, it is reasonable to assume a power 
law, 

^ oc ES. (16) 

dE e ± e± v ; 

The normalization can then be extracted by a simple integration above the relevant energy 
threshold, conventionally taken to be ~ 1 GeV. Information on the spectral index a can be 
inferred from actual multi-wavelenght pulsar observations. Secondary synchrotron radiation 
at radio frequencies indicate, f or th e e ± population of e.g. pulsar wind nebulae, spectral 
indexes ranging from 1 to 1.6 195[. Larger values of a, or a steepening of the spectrum 
at larger energies, are implied from the detected X-ray emission from the same objects. In 
addition, e.g. in the HR model the spectral index of gamma-ray from pulsars should be 
close to the spectral index of the injected e ± . EGRET observations of galactic pulsars point 
;o fa irly 



196 



lard gamma-ray spectra in the 0.1 to 10 GeV range, loosely between 1.4 and 2.2 
1971 ] . For all these reasons, we shall consider here the rather liberal range of values 
1.4 < a < 2.4. 

In addition to the electron-positron spectrum, a potentially important impact on the 
measured e ± flux on Earth from pulsar physics stems from the finite time during which 
the energy output estimated above is delivered in the inter-stellar medium. As shown e.g. 



in fig. 4 of Ref. 23|, different decay time scales, or different models for the injection 
time profile, impact the spectrum of the produced in pulsars. Specifically, a tail at 
high energy extends in the e ± spectrum beyond the cutoff energy set by inverse Compton 
and synchrotron losses. There is no reason to believe that there exist one unique injection 
spectrum for galactic pulsar emission of electrons and positrons. Including this effect is 
ultimately important (and technically very easy, once the form of the light curve is assumed), 
but it would only add an extra parameter (per pulsar) to the present discussion. As for the 
injection spectrum cutoff, adding more parameters would only make it easier to demonstrate 
the point we want to make: existing known pulsars provide a natural explanation to the 
positron fraction and to the electron-positron differential spectrum data. Therefore, we 
shall hereafter assume a bursting e ± injection, i.e. that the duration of the emission is much 



shorter than the travel time from the source 35[, which, incidentally, is an appropriate 
approximation for almost all pulsars we consider here. 

After production, the e ± emitted from a galactic pulsar diffuse and loose energy prior 
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to arrival at the Earth's atmosphere. These processes have been described in the standard 



diffusion approximation (i.e. neglecting convection) by a familiar diffusion equation [16], 
which, for local sources (in a leaky box type scenario, where the height of the diffusion 
region is larger than the sphere where the contributing sources lies) can be assumed to have 
spherical symmetry, and reduces to the form 2^|: 

df D( 7 ) d 2 8f d{P(i) f) „ 
at r A or or 07 

In the equation above, /(r, i, 7), with 7 = E e ±/m e , is the energy distribution function of 
particles at instant t and distance r from the source; 

P(l) = Po + Pil + P2I 2 (18) 

indicates the energy loss term (which for our purposes can be very accurately approximated 
by the last term only, describing the dominant inverse Compton and Synchrotron losses, 
P(7) ~ P27 2 at high energy); D{pf) oc 7 5 is the energy-dependent diffusion coefficient, which 
is here assumed to be independent of r (i.e. we assume a homogeneously diffusive medium; 



for a discussion of the role of a non- homogeneous medium, see App. C in Ref. 38[). 

As mentioned above, the dominant e p m energy losses at energies at or above E e ± > 10 
GeV are (i) inverse Compton scattering off photon backgrounds such as the cosmic microwave 
background (CMB) and photons at optical and IR frequencies, and (ii) synchrotron losses 



in the ISM magnetic fields. Following Ref. 



231 ] . we estimate: 



P2 = 5.2 x 10- 20 — ^° b-\ (19) 
1 ev/cm 

where 

w = w B + wcmb + w opt (20) 
and (again following the estimates of Ref. 23]) where u>cmb — 0.25 eV/cm 3 and w op t — 



0.5 eV/cm 3 indicate the CMB and optical-IR radiation energy density, while wb 



0.6 eV/cm 3 (for an average B = 5/iG) is the energy density associated to galactic mag- 



netic fields. Altogether, we rather conservatively assume P2 — 5.2 x 10 20 s 



Under the assumption o: 
readily calculated to be 



a power law injection spectrum, the solution to eq. (ITTj) can be 



^^) = §^(l-p 2 t7r 2 (^Vexp(-(-) ], (21) 
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with 7 < 7 cut = (p2t) 1 , and / = otherwise, and where 



r »"f |m l W ' (22) 

Assuming a pulsar burst-like injection at t — T, where T is the age of the pulsar, the local 
spectrum of e from pulsars is entirely fixed once the distance Dist = r to the pulsar, its 
age and the normalization to the injected spectrum are given. THis approximation is well 
justified for mature pulsars, where T is much larger than the trapping time of before 
they diffuse in the ISM, given the steep breaking index of dipolar emission as a function of 
time. The local e flux will then simply be: 

J = cxf{Dist,T 1 y)/(4ir). (23) 

The first important feature of eq. ( 12~TI) above is that the pulsar spectra have a sharp cutoff 
(possibly smoothed out by deviations from the power-law behavior in the injected at high 
energy, expected in the very high energy regime, and by non-burst-like injections), and a 
normalization that depends only on the (energy dependent, unless 5 = 0) factor s = r/r^. 

As clear from the discussion above, the pulsar distance is therefore, clearly, a parameter of 
crucial importance in the estimate of the flux of the pulsar-injected electrons and positrons. 
For consistency, we refer to all pulsar distances as listed in the ATNF catalogue. Since, 
later on, we will further discuss the case of Geminga, and given that the estimate of the e 
flux sensibly depends upon distance, we comment further on the value of 160 pc listed in 
the ATNF catalogue. This value refers to the result of optical measuremen ts of the annual 
parallax of Geminga with data from Hubble Space Telescope observations 198], yielding a 
distance of 157^4 pc. T he range is compatible with the more recent and more conservative 



result reported in |l99j . 250^2° P c - The large error-bars in the last measurement should 
serve as a warning to the reader as to the significant uncertainty in the measurement of 
pulsar distances, even in the case of well-known and studied objects like Geminga. Distances 
inferred from dispersion measures rely on the measurement of the mean electron number 
density at di fferen t places in the Galaxy, and are believed to be affected by errors as large 



factor 2 [1861 ] . 



In addition to uncertainties in the current distance of pulsars, an important role in the 
estimate of the cosmic-ray electrons and positrons injected by mature pulsars is played 
by the pulsars' velocities. Pulsars have an average velocity at birth of around 250-300 
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TABLE I: The three diffusion setups we consider in the present analysis. D\ indicates the value 
of the diffusion coefficient at E e ± = 1 GeV, in units of cm 2 /s. 

Scenario D\ 5 

MAX 1.8 x 10 27 0.55 

MED 3.4 x 10 27 0.7 

MIN 2.3 x 10 28 0.46 



km/s 200], and a distribution of observed velocities today consistent with a mawellian 
distribution with a large width, of around 200 km/s. In addition, several outliers have been 
observed, with velocities even in excess of 1000 km/s. This means, for instance, that a 
pulsar with the current distance of Geminga, and with an age close to that estimated by 
its timing parameters, could have been almost anywhere in the sky at birth, and hence 
when it inject ed in the ISM the high energy e ± population we are interested in here. For 



instance, Ref. |198| points out that if Geminga has a large enough transverse velocity (> 700 
km/s), it could have traveled, since its birth, 250 pc and have been potentially located, at 
birth, at coordinates diametrically opposite to those observed today. The issue of pulsar 
proper velocities is clearly a very relevant one in the game of determining if a given pulsar 
is the origin of the observed lepton anomalies. This issue is particularly acute for nearby 
and relatively old pulsars, while is less severe for more distant pulsars. In this respect, the 
analysis we present below implicitly assumes, for nearby pulsars, a small transverse velocity. 
We show quantitatively the effect of pulsars' velocities on the pulsar age-versus-distance 
parameter space in fig. (TUJ 

The last ingredient before starting our estimates is to define the diffusion coefficient. 
Since we only care about E e ± > 10 GeV, we neglect not only the effects of solar modulation, 
but also the energy- independent part of -D(7) that kicks in at low energies, and assume 
D{^) oc 7 5 , with the combination of normalizations and values of 5 specified in tab. (H 
for the three diffusion setups MIN, MED and MAX. These three scenarios are such that 
simulations of various cosmic ray species abundan ces w ith those setups reproduce available 



cosmic-ray data, according to the analysis of Ref. [201 ] . 

Notice that the analytical solution of Eq. (]2ip assumes spherical symmetry, an assumption 
which can clearly break down for the geometry of the diffusive region of the Galaxy, ordinarily 
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assumed to be cylindrical, with a diffusive region of height L. Specifically, the three diffusion 
models of Tab. [J correspond to half-thiknesses of the diffusion zone L = 1, 4 and 15 kpc. 
Imposing free escape boundary conditions, namely that the cosmic ray density equals zero at 
the surfaces of the diffusive region (i.e. on the planes at z — ±L) can b e en forced exactly by 



resorting to a set of image charges similar to the one outlined in Ref. [2021 ] . Specifically, we 
need to add an infinite series of image-pulsars identical to the one under consideration (that 
injects at a time T, with a given injection spectrum etc.) at z n = 2Ln + (— l) n £cb where 
Zq is the distance of the original pulsar fro m the galactic plane, and adding all image-pulsars 



(for — oo < n < +oo) with a factor (— l) n [2021 ] . Assuming z ~ 0, this amounts to simply 



account for the suppression factor cyl(£) given by 

cyl(£) = Yl (-!)" ex P K v ) > (24) 

n=— oo 

where £ defined as 

In the general case zo ^ 0, we simply have 

+oo 

cyl(£, zo) = ex P ~ 2 (- 1 )^f ) • (26) 

n=— oo 

The Reader should bear in ming that the method outlined above neglects the radial boundary 
condition usually considered in diffusion models. This is well justified if the sources under 
cons ideration are closer to the observer than to the radial boundary of the diffusive region 



202| . This is a good approximation with the galactic sources under consideration here (see 



e.g. fig. 10.5 of [186J). Interestingly, cyl(£) doesn't depend on the pulsar distance, but only 
on the ratio of the thickness of the diffusion zone over the effective diffusion radius rdif. 
Quantitatively, this statement holds as long as the pulsar distance is small compared to the 
neglected radial boundary condition. 

Fig. [3] illustrates the effect of deviations from the spherically symmetric solution as a 
function of £, as well as for selected pairs of pulsar ages and e energies, as a function of L 
(we employ the values of the MED parameters for this plot, but let L vary). As a rule, the 
breakdown of spherical symmetry must be taken into account when £ < 0.5, and a source 
becomes essentially irrelevant if £ < 0.2. Effects of order 10% are expected for diffusion 
regions with a thickness between 1 and 2 kpc, and for the high-energy end (E e ± ~ 10 GeV) 
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FIG. 3: The effect of electron-positron escape from the Galaxy (defined by the suppression fraction 
cyl(£)), as a function of the thickness of the galactic diffusion region L, for selected pairs of values 
for the injection time and energy. 

of the spectrum of young pulsars, or for the low energy end (E e ± ~ 10 GeV) of mature 
pulsars. For the following discussion, escape of pulsar- injected e ± is therefore typically not 
relevant. 

To summarize, once a theoretical model is chosen, the contribution of pulsars to the 
flux depends on three observational inputs: the pulsar age (as determined from timing 
measurements), its distance and its spin-down power (which can be connected, assuming a 
particular model such as one of those outlined here above, to the total e energy output 
injected in the ISM). To compute the resulting flux at Earth, one further needs to specify 
a diffusion setup: we restrict this choice to the three possibilities of tab. [B We are now 
ready to calculate the contribution of some of the most likely sources of primary energetic 
positrons (and electrons): nearby bright pulsars. 
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TABLE II: Data for a few selected nearby pulsars and SNR's. E OVLt is the energy output in e 
pairs in units of 10 48 erg (for the ST model column we assumed f e ± = 3%). The energy output 
for the SNR Loop I and Cygnus Loop are not estimated within the ST model, but via estimates 
of the total SNR output. The f e ± column indicates the output fraction used to compute the 



fluxes shown in fig. 0] and [5] assuming the ST model. 



Name 


Distance [kpc] 


Age [yr] 


E [ergs/s] 


E ut [ST] 


Bout [CCY] 


-Eout [HR] 


Eout [ZC] 




a 


Gcminga [J0633+1746] 


0.16 


3.42 x 10 5 


3.2 x 10 34 


0.360 


0.344 


0.013 


0.053 


0.005 


0.70 


Monogem [B0656+14] 


0.29 


1.11 x 10 5 


3.8 x 10 34 


0.044 


0.133 


0.006 


0.020 


0.020 


0.70 


Vela [B0833-45] 


0.29 


1.13 x 10 4 


6.9 x 10 36 


0.084 


0.456 


0.006 


0.372 


0.0015 


0.14 


B0355+54 


1.10 


5.64 x 10 5 


4.5 x 10 34 


1.366 


0.677 


0.022 


0.121 


0.2 


0.61 


Loop I [SNR] 
Cygnus Loop [SNR] 


0.17 
0.44 


2 x 10 5 
2 x 10 4 




0.3 
0.03 








0.006 
0.01 





i i i i i i i i i i i i i i i i i r 




Energy [GeV] 

FIG. 4: The spectrum of selected nearby pulsars and SNR's (for the parameters employed to 
calculate the fluxes see tab. |TTJ) . We assume an e ± injection spectral index a = 2, and a median 
diffusion setup (MED). The dotted lines correspond to injection spectra featuring an exponential 
cutoff at E e ± = 10 TeV. 
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- 0C=2, MED diffusion setup, ST model 




Energy [GeV] 

FIG. 5: The positron fraction for the same sources as in fig. [Hand tab. HH 
III. CONTRIBUTIONS FROM SELECTED NEARBY PULSARS 

In this section we apply the setup outlined above to a few objects that, naively, can 
contribute to the e ± flux needed to account for the PAMELA and ATIC data. Specifically, 
we consider the Geminga [J0633+1746], the Monogem [B06546+14] and the Vela [B0833- 
45] pulsars: all these pulsars lie within 300 pc, and although Vela might be too young to 
contribute significantly to the flux in the energy range of interest, and the produced e ± 
might not have diffused in the ISM yet, all these pulsars have a very large spin-down energy 
loss rate. In addition, we consider the radio pulsar PSR B0355+54, a middle aged (T ~ 0.56 
Myr) pulsar with a fairly large spin-down energy loss rate and a relatively small distance. 
We are lead to consider this particular pulsar to explain the ATIC "bump" feature, i.e. the 
reported excess over the expected e ± flux. Naively, a pulsar that contributes to the bump, 
assuming a relatively hard injection spectrum, should have an age 

Tatic ~ m e / {p 2 x £ bump ) « 5 - 6 x 10 5 yr (27) 
for E hump ~ 500 - 600 GeV. (28) 
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In addition, the candidate pulsar should be distant enough for the resulting spectrum to be 
peaked enough around E^ nmp to avoid excessive contributions at lower energies, and bright 
enough to produce the right level of e ± . In short, PSR B0355+54 appears to be a natural 
candidate fulfilling all of these requirements, which is why we consider it here (notice that 
in what follows we consider the possibility of multiple pulsars explaining the ATIC feature 
as well). Finall y, we also include the two SNR's Loop I and Cygnus Loop, following the 



analysis of Ref. 



203]. As mentioned in the introduction, radioactive decay, or even a faint 



and yet undetected pulsar associated to these SNR could inject enough positrons for the two 
SNR's to contribute significantly to the positron fraction as well. 

We outline the set of pulsars for which we shall estimate the positron fraction and the total 

spectrum in tab. [Ill The table indicates the name, distance, age and the spin-down power 
of the various pulsars. The following four columns specify the estimated total e ± energy 
output according to the four models described in sec. HH with our standard assumptions for 
the relevant model parameters. Finally, the f e ± column indicates the actual value for the 
efficiency factor in the output within the ST model. In all cases we adopt values around 
the estimate quoted in Ref. 35J and in Ref. 38|. In one case, for PSR B0355+54 we adopt 
a factor 6-7 larger efficiency output. This should be contrasted with the variance in the 
model estimates for the total output of pulsars, but it certainly indicates an unusually 
large efficiency for this single pulsar to produce enough e ± to account for the ATIC bump. 

Fig. |4] shows the resulting differential energy spectrum (times energy to the third 
power) for the six sources listed in tab. [Til as we h as the ATIC Q and the HESS data lo) ]. 
The following fig. [5] shows the estimated positron fraction, with the PAMELA data jl|] . We 
assume a median spectral index a = 2 for all sources, and employ the MED diffusion setup 
and the ST output model, with the values of f e ± specified in tab. [II] For illustrative 
purposes, we also show, with dotted thin lines the effect of a cutoff at E cnt = 10 TeV for the 
Monogem and Vela pulsars. 

The figures show that both the nearby Geminga pulsar and the Loop I SNR, both at 
distances around 160-170 pc, and both with ages T ~ 2 — 3 x 10 5 , can very naturally be 
the dominant positron sources to explain the PAMELA data, with an output in less 
than a percent of the pulsar rotational energy, in line with, or less than, previous estimates 
of that efficiency factor. The resulting contribution to the total e ± at larger energies is 
subdominant. Young, powerful and relatively nearby sources like Vela and the Cygnus 



25 



Loop likely contribute to the highest energy part of the total differential spectrum. 
Fig. H] seems even to suggest that the highest energy H.E.S.S. bin could even actually reflect 
the contribution from these bright local sources. The contribution of these sources to the 
positron fraction kicks in above the TeV scale, and is therefore irrelevant for the positrons 
detected by PAMELA. The Monogem SNR pulsar might plausibly contribute to the higher 
energy bins reported by PAMELA, and could be dominating the e^ 1 flux between 1 and 3 
TeV, even with a efficiency factor of f e ± ~ 1.5%. Finally, the ATIC bump seems to fit 
very naturally with the expected e flux from the PSR B0355+54, although this requires 
a very efficient, but still possible in principle, output of (/ e ± ~ 20%). The positrons 
produced by PSR B0355+54 would give a sizable and steeply rising contribution to the 
positron fraction between 200 GeV and 600 GeV, a range that might eventually be explored 
by PAMELA Q. 

In short, in this section we showed that selected, known nearby e^ 1 sources, with a reason- 
able e ± energy output, very naturally account for both the PAMELA and the ATIC data. 
None of these sources appears to be suitable to explain both experimental results at once. 
In the next section we address the question of whether one single e source (combing the 
entire pulsar parameter space) can explain at once the PAMELA and the ATIC data. 

IV. MINIMALISM: SINGLE SOURCE 

In fig. [6] and [7J we explore the possibility that the ATIC and PAMELA data originate 
from one single source. We fix the injection spectrum of at the source to a =1.4, 1.6, 1.8 
(fig. [6]) and to a =2.0, 2.2, 2,4 (fig. [7]). For each value of a, in the top, middle and bottom 
panels we adopt the MAX, MED and MIN diffusion setups, respectively. On each panel, 
we normalize the source e ± energy output to get the best fit (i.e. the lowest \ 2 P er degree 
of freedom (d.o.f.)) to both the ATIC and the PAMELA data. The regions shaded in grey 
have fluxes, in at least one of the H.E.S.S. data bins, in excess of the measurement by 2-cr, 
and we therefore consider these regions excluded. The black contours indicate a global \ 2 
per degree of freedom less or equal than 3. The best fit regions are therefore inside the black 
contours. Notice that we verified that including other data sets on the high-energy flux of 
electrons and positrons (e.g. the recent PPB-BETS measurement j^j]) does not visibly affect 
the contours corresponding to the ATIC data alone: this simply depends on the much better 
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FIG. 6: Likelihood contours for a single pulsar-like source normalized to give the best joint fit 
to the PAMELA and to the ATIC data, in the age versus distance plane. The best fit regions 
(x 2 /d.o.f. < 3) are below the black curves. The regions within the red and green lines have a 
X 2 /d.o.f. < 1 for the fit to the ATIC and to the PAMELA data alone, respectively. The grey 
region is excluded by the HESS data [lj|. The blue contours indicate the inferred energy 
output, in units of 10 48 erg. The black circles indicate the age and distance of known pulsars in 
the ATNF catalogue. In the three upper panels we assume the MAX diffusion setup, in the central 
ones the MED setup and in the three lower panels the MIN setup. The three panels to the left 
have e ± injection spectral index a = 1.4, those in the middle 1.6 and those to the right 1.8. 



statistics that the ATIC measurements have with respect to previous data. 

The regions inside the green and the red contours have a \ 2 P er degree of freedom less or 
equal than one for the individual ATIC and PAMELA data-sets, respectively. This means 
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FIG. 7: As in fig. [6j for softer injection spectral indexes (2, 2.2 and 2.4, from left to right). 

that, by normalizing the source to the best fit to both data sets, the regions giving the best 
fits to the individual ATIC and PAMELA data are those within the green and red contours. 

Finally, we show contours of constant energy output, which we indicate as W = 
E out /10 48 erg. Roughly, a reasonable range for W is between 0.1 and 10. We also indicate 
with black circles the actual location of the ATNF pulsars in the age versus distance plane. 
We, again, stress that those points locating pulsars in the age-distance parameter space are 
subject to (i) uncertainties in the determination of the current distance of the pulsars, and 
on the determination of the age (what is shown is the characteristic age estimated from time 
measurements), and (ii) a systematic uncertainty from the individual pulsar velocities (this 
effect is particularly relevant for nearby and mature pulsars, hence for those points in the 
lower-right corner of each panel). 

We also remark that the contours we show in fig. [6] and [7] would be naturally affected 
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by changing the astrophysical background model for the diffuse cosmic-ray electron and 
positron component. Specifically, we concentrate on what would happen if one changed the 
normalization and the spectral index of the diffuse galactic high-energy electron-positron 
spectrum. We verified that a shift in the local normalization almost only affects the contours 
by shifting the contours of iso-level for the variable W, but does not affect significantly the 
shape of the regions preferred by an explanation to the PAMELA and ATIC anomalies. A 
tilt in the diffuse background does affect the shape of the contours. Remarkably, though, this 
effect amounts to shifting the contours towards harder or softer values for a, that compensate 
the variation in the background spectrum. In essence, we find that changing the background 
spectrum corresponds to variations in the location and shape of the contours favored by the 
lepton anomalies comparable to what shown, from left to right, in fig. [6]and[7J 

With a hard spectrum, fig. [6], clearly a single source interpretation is disfavored: the 
ATIC data strongly prefer a pulsar with an age around 0.4 to 0.6 Myr, while PAMELA 
favors an older source with an age in the 3-6 Myr range 3 . As the spectrum gets softer, the 
PAMELA data start to be compatible with a broader combination of ages and distances. 
For a = 1.8, the Geminga pulsar appears to be roughly at the intersection of the PAMELA 
and ATIC xVd.o.f. = 1 regions. 

Moving to softer injection spectra (fig. [7J, the ATIC data favor a very distinct region in 
the age-distance plane, again with ages around 0.4 to 0.6 Myr, and distance ranges which 
depend on the diffusion setup and on the injection spectral index. PAMELA points toward a 
relatively wide range of very correlated values for age and distance. Several pulsar candidates 
fall within both the ATIC and the PAMELA xV d -°- f • = 1 candidates, but again, the overall 
information we extract is that the two regions are fairly distant. The ATIC and PAMELA 
data seem to most naturally require two distinct sources. 

As a last comment, looking at the blue lines with contours of output power, the regions 
favored by ATIC tend to require very large e energy outputs, especially with softer injection 
spectra, pointing to the fact that the pulsar scenario is compatible, but does not favor, for 
natural values of the parameters, a peaked spectrum such as what reported by ATIC. 

Given that a single source interpretation gives only a partly satisfactory fit to the data, we 

3 Notice that the precise value of these age ranges depends critically on the assumed e ± quadratic energy 
loss coefficient P2- 
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now move on to explore the region of pulsar parameter space which are, separately, favored 
by the PAMELA (sec. and by the ATIC (sec. ED> data. 

V. PAMELA AND OCCAM 

The four panels of fig. [8] outline, for various pairs of injection spectral indexes and diffusion 
setups, the region of the pulsar age versus distance parameter space favored by PAMELA 
and compatible with the ATIC and with the H.E.S.S. data. 

The two panels at the top assume a MED diffusion setup with e^ 1 injection spectral index 
a = 1.6 and 2.4, while the two panels at the bottom employ a MAX diffusion setup, with 
a = 1.8 and 2.2. As in the figures above, the grey regions are in conflict with the H.E.S.S. 
data, and the blue lines indicate the source energy output in e ± , in units of 10 48 erg. In 
addition, we shade in orange the regions disfavored by the ATIC measurements, and label 
some of the pulsars which might be relevant to explain the PAMELA signal. The dashed 
contour outlines the region inside which x 2 /d.o.f < 1, while inside the black solid contours 
X 2 /d.o.f < 0.5. 

With a hard a ~ 1.4... 1.6 spectrum, the PAMELA data point to nearby pulsars (distance 
less than 1 kpc), with an age (correlated to the distance) which the H.E.S.S. data constrain 
to be larger than Ri4x 10 5 yr, and that the PAMELA data force to be less than « 4 x 10 6 yr. 
Consistency with the ATIC data, as well as the values for the energy output, point towards 
a multiple source contribution as the most likely situation. The first four entries of tab. II III 
list a few plausible pulsars that can contribute to the PAMELA signal, assuming a MAX 
diffusion setup. The list includes PSR J1918+1541, B0450+55, B0834+06 and B1845-19. 

Softer spectra generically worsen the quality of the fit to the PAMELA data (of course, 
not all pulsars will have the same spectral index, but softer spectra are generically unlikely 
to add up to the positrons measured by PAMELA with the correct spectral shape). With 
the MED diffusion setup, the Monogem pulsar, with a reasonable e ± energy output close 
to ~ 10 48 erg, gives an acceptable fit to PAMELA. With a MAX diffusion setup, it instead 
seems more plausible that a series of pulsars with lifetimes of a few Myr and distances 
between 0.7 and 1 kpc are the most relevant contributors to the PAMELA signal. 

In any case, we showed in this section that no matter which injection spectrum is as- 
sumed, and which diffusion setup is employed, known pulsars naturally predict a positron 



30 



Best Fit to Pamela data, "MED", a=1.6 



Best Fit to Pamela data, "MED", a=2.4 



1000 



100 



• 1 11 


4 , ,,, t . 

Excl. by 
*. 

ATJ£^ 


• • • * « 


•T 1 1 
• * 

• 


Excluded by HESS (2a) ^ 

; ^065.3 W=10/ 
HB21 / 

/ W=l 


/.B/)jkif9 * 
B0458T+55'' 1,' V * 

Ji9.f8+js4i«/* | «\ 

. /TC>n<*t-i<; l. 


Cy^nus Loop 

• 7 / ' 

/ / / 
/ / / 
/ / / / 

/Vote B065>$14/ 

y w=o.i • / / . 


/ /\/ 


-3125 1 
■ * 1 • 




/ >/ / 
/ / / / 

- ' / Y 

/ / / /^voop I 

' / /* \ • 

/ / \ Geminga 

/ / / 1 

SJ,/a , 1 m 


A 

i 

/J1856-3754 
• 





10 



10 



10 



10' 



Age [yr] 

Best Fit to Pamela data, "MAX", ot=1.8 



1000 




Age [yr] 



1000 



100 



• 1 1 1 * "V 

Excl. by 


. . 1 1.1^11 


• # • 

• •/ '••* »' 
•v / A *' 


* ATIC 




Excluded by HESS (2c) \ y 




- • / 




G65.3 / / • • 


HB21 / / 

/ / 


• • - 

• • 


// 


• 


Cygnus Loop ^ ^/ 


• 


/ / J0T20 


-3125 


' \ 
Vela ."B0656<14 \ 


• • a 




/ 




/ \ w 


= 10 


, / noop I 
^ / / • \ 




/ \ • 

* y \ Geminga 




/ / \ 


J1856-3754 


/ W/=0.1 \\|=1 


• 



10 



10' 



10 



10' 



Age [yr] 

Best Fit to Pamela data, "MAX", a=2.2 



1000 



Excluded by HESS (2a) 

G65.3 
m 

HB21 

Cygnus Loop 



• • * 



Excluded by • *\ • ^ 
ATIC B045,rf+55l • r 

•/ • 




Age [yr] 



FIG. 8: Regions providing the best fit, with a single source, to the PAMELA data The grey 
region is ruled out by HESS, and the orange region is in conflict with the ATIC data. As in fig. [6J 
we indicate (and sometimes label) the age and distance of known pulsars, as well as curves of 
constant inferred e 1 * 1 energy output, in units of 10 48 erg. The region of parameter space inside the 
solid black lines has a x 2 /d.o.f. < 0.5, while inside the black dashed region % 2 /d.o.f. < 1. 

fraction close to what is reported by PAMELA, with very reasonable electron-positron en- 
ergy outputs. 
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VI. ATIC AND OCCAM 



Fig. [9] shows the portion of the pulsar age-distance parameter space favored by the ATIC 
data. The notation is as in the previous figure [HI but now the regions shaded in cyan 
are disfavored by PAMELA. Independently of the injection spectrum, the inverse Compton 
plus Synchrotron cutoff, as mentioned above, dictates an age between 0.4 and 0.6 Myr. In 
addition, while distance matters less for a hard spectral index (as long as the pulsar e ± 
output is large enough), for softer spectral indexes the ATIC data prefer distances of the 
order of 1 kpc. 

An important message we read off of fig. M is that unless the injection spectrum is very 
hard (left panels), existing pulsars must be extremely powerful to account for the ATIC 
data: in the right panels the pulsars that lie within the x 2 /d.o.f < 1 regions typically need 
normalizations W ~ 100, i.e. an e ± energy output of E out ~ 10 50 erg, which, while possible 
in principle in some of the models we considered above (see fig. [1]) seems unlikely for the 
standard energy budget of supernova explosions. 

In summary, while the PAMELA data can be regarded as a natural posi-diction of the 
pulsar scenario, the ATIC data find a possible explanation in this scenario as well, although 
energy output considerations indicate that the signal reported by ATIC would not be the 
most obvious to be expected from known pulsars. 



VII. ONE SINGLE, NEARBY SOURCE VERSUS A COMBINATION OF PUL- 
SARS 



In this section we consider the coherent superposition of a few pulsars that could con- 
tribute to the ATIC signal (namely, PSR B0919+06, B0355+54, B1055+54, J1849-0317 and 
B1742-30, see fig. [9]), together with the objects, listed above, that might contribute to the 
PAMELA signal. We give the age, distance and the estimated e energy output for these 
pulsars in tab. IIHI Guided by the analysis we carried out in the two previous sections, we 
want to contrast, here, a scenario where the PAMELA and ATIC data are accounted for 



by a single bright source, namely Geminga 35j, |36[, to a scenario where a combination of 



pulsars coherently contributes to the observed e ± fluxes (tab. IIIII) . 

We adopt the MAX diffusion setup, and a spectral index a = 1.85 for all pulsars, including 
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FIG. 9: Regions providing the best fit, with a single source, to the ATIC "bump" feature [2] (we 
define the bump feature as all ATIC data with E e ± > 300 GeV). The grey region is ruled out by 
HESS, and the cyan region is in conflict with the PAMELA data. As in fig. we indicate (and 
sometimes label) the age and distance of known pulsars, as well as curves of constant inferred 
energy output, in units of 10 48 erg. The region of parameter space inside the solid black lines has 
a x 2 /d.o.f. < 0.5, while inside the black dashed region x 2 /d.o.f. < 1. 
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TABLE III: Possible combinations of multiple pulsars contributing to explain the PAMELA and 
the ATIC data. P/A refers to whether the pulsar dominantly contributes to the PAMELA or to 
the ATIC signal. E ou t is the energy output in pairs in units of 10 48 erg. 



Name 


P/A 


Distance [kpc] 


Age [yr] 


Km [ST] 


E ont [CCY] 


Km [HR] 


E ut [ZC] 


J1918+1541 


P 


0.68 


2.31 x 10 6 


0.99 


0.33 


0.023 


0.022 


B0450+55 


P 


0.79 


2.28 x 10 6 


1.16 


0.37 


0.025 


0.025 


B0834+06 


P 


0.72 


2.97 x 10 6 


0.11 


0.07 


0.011 


0.001 


B1845-19 


P 


0.95 


2.93 x 10 6 


0.01 


0.015 


0.005 


0.0002 


B0919+06 


A 


1.20 


4.97 x 10 5 


0.158 


0.178 


0.010 


0.016 


B0355+54 


A 


1.10 


5.64 x 10 5 


1.366 


0.677 


0.022 


0.121 


B1055-52 


A 


1.53 


5.35 x 10 5 


0.82 


0.49 


0.017 


0.075 


J1849-0317 


A 


1.90 


4.81 x 10 5 


0.06 


0.10 


0.007 


0.007 


B 1742-30 


A 


2.08 


5.46 x 10 5 


0.24 


0.22 


0.012 


0.022 



Geminga. In addition, we use the ST model with the reference positron-electron output 
efficiency factor f e ± = 3% for the combination of pulsars, and a lower efficiency facto r for 
the single source Geminga, f e ± = 0.45%. Also, following the estimate of Ref. 203], we 
assume an exponential cutoff for the Geminga injection spectrum, at an energy of 0.67 TeV. 
It is of course a constraining assumption to take the same spectral index and the same 
positron-electron output efficiency for all pulsars in the combination we adopt. There is no 
reason to believe that every pulsar has the same value for these parameters. However, if 
we indeed find satisfactory agreement with data, this will reassure us and make our point 
stronger; had we not found good agreement, we could have massaged the individual pulsar 
parameters and gotten much better agreement. 

We show in fig. [10] the location on the age versus distance plane of the pulsars reported in 
tab. [Jill which we use as combinations favored by observational data. Again, we point out 
that the ATIC data appear to be compatible with an origin from sources located at a distance 
of 1-2 kpc, and in a narrow range of ages, 0.4-0.6 Myr. PAMELA points towards closer (0.7-1 
kpc) and older (2-4 Myr) pulsars instead. We also show t he po sition of the Geminga pulsar. 
For the standard background we use the best GALPROP [2 04 ] self-consistent predictions for 
the primary and secondary electrons plus the secondary positrons, and we add a cutoff in the 
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FIG. 10: The location, in the age versus distance plane, of the pulsars of tab. IIIIl The vertical 
arrows indicate the potential, maximal effect of the pulsars' proper velocity, for a velocity of 300 
km/s (solid lines) and of 1000 km/s (dashed lines). The horizontal arrows pointing outwards from 
the Geminga pulsar point indicate a range of a factor 2 uncertainty in the determination of the 
actual pulsar age from the characteristic age obtained from timing measurement (dotted lines) and 
the effect of the finite time it takes for the injected by the pulsar to diffuse in the ISM, assumed 
to be of order 10 5 yr (solid line). 



primary electron spectrum around 1 TeV, compatibly with the H.E.S.S. data of Ref. [Lj 

In fig. [TU]we wish to draw again the attention of the reader to the uncertainties that affect 
the determination of the distance and age of pulsars when they actually inject high-energy 
cosmic-ray electrons and positrons in the ISM. In particular, the vertical arrows indicate the 
potential, maximal effect of the pulsars' proper velocity, for a velocity of 300 km/s (solid 
lines) and of 100 km/s (dashed lines), values close to the average and maximal range of 
pulsar kicks 2001 ] . The horizontal arrows pointing outwards from the Geminga pulsar point 
visually indicate, instead, a range of a factor 2 uncertainty in the determination of the actual 
pulsar age from the characteristic age obtained from timing measurement (dotted lines) and 
the effect of the finite time it takes for the e ± injected by the pulsar to diffuse in the ISM, 
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FIG. 11: The positron ratio (left) and the e ± flux (right) for the Geminga pulsar (with f e ± = 0.45%) 
and for the combination of multiple pulsars specified in tab. Ill II (with f e ± = 3%). We assume the 
MAX diffusion setup and an injection spectral index a = 1.85. 

assumed to be of order 10 5 yr (solid line). 

The resulting positron fraction and overall differential energy spectrum for the com- 
bination of pulsars and for Geminga are shown in fig. [TTJ As far as the PAMELA data 
are concerned, we obtain in both cases excellent fits. Again, had we chosen a non-trivial 
distribution of spectral indexes in the pulsar combination, the agreement could have been 
arbitrarily better. The positron fraction in the two cases (single versus combination) has 
vastly different features: The "ATIC pulsars" induce a sharp rise in the positron fraction 
after the abrupt cutoff of the "PAMELA" pulsars. The expected positron fraction therefore 
has a two-bumps structure. Of course, there is no reason not to assume that more pulsars 
would contribute at intermediate energies (E e ± ~100...200 GeV) and somewhat fill in the 
gap. In the case of Geminga, the prediction is a sharp drop around one TeV, and a consistent 
growing trend with energy below a TeV. 

The right panel shows the overall e spectrum in the two scenarios. Clearly, a single 
source does a rather poor job at fitting the detailed features in the ATIC data, while a 
combination of pulsars traces almost perfectly the feature reported by that experiment in 
the 600 GeV range. Notice that a combination of pulsars also shows a prominent drop-off 
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feature at lower energies, which could be detected with higher statistics, e.g. by the Large 
Area Telescope (LAT) onboard Fermi (see sec. IIXI) . 

In summary, the PAMELA data are perfectly compatible with both a single source and 
with a combination of coherently contributing pulsars. The ATIC data favor a scenario with 
more than one contributing source, and the resulting combination produces an interesting 
and very clean signature in the positron fraction. 

In this section we considered only a sub-section of the numerous pulsars in the ATNF 
catalogue. It is important and crucial, however, to understand whether our predictions are 
consistent with adding up all other objects in the catalogue, and where and how the more 
than 1,000 other pulsars contribute to the production of energetic e ± . We address this items 
in the next section. 



VIII. SANITY CHECK: ADDING UP ALL ATNF SOURCES 



We study in this section the overall contribution of all g < 1 ("gamma-ray") relatively 
distant pulsars (distances larger than 1 kpc) to the positron fraction and to the e differential 
flux. Overall, we include 248 objects, of the 266 total gamma-ray pulsars in the ATNF 
catalogue. The remaining nearby pulsars (distances less than 1 kpc) which are not included 
here have either been discussed in sec. IHII or in sec. |V] IVII and [VH[ or are too faint to give 
any significant contribution. 

Scope of this section is to ascertain that within the models outlined in sec. [Til for the 
pulsar e ± energy output, existing pulsars do not over-produce energetic e ± . This exercise is 
therefore a "sanity check" to our approach. On the other hand, we will be able to appreciate 
what the contribution of existing, well-known pulsars is, for various diffusion models and 
average injection spectral index, and to decompose it in terms of distance and age. 

While our approach is closer to observations than theoretical extrapolations based on 



assumed spatial distributions and pulsar birth rates |35|, |37| , an important point to keep in 
mind is the issue of incompleteness: as for all astrophysical catalogues, not all galactic pulsars 
have been detected yet (and are therefore not included in the ATNF catalogue) within the 
age-distance domain we consider here. One aspect of the catalogue incompleteness is due 
to the simple fact that some fraction of the pulsars that could contribute to the local flux 
of are not bright enough to have been detected so far. A second aspect is that the 
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FIG. 12: The contribution from all gamma-ray pulsars (g < 1) in the ATNF catalogue located at 
a distance greater or equal than 1 kpc. The left panels show the positron fraction, and the right 
panels the flux from the pulsars. In the torj^panels we study the variation in the contribution 
from relatively distant sources from the model with which we compute the output. In the 



pulsars radio emission is beamed, and depending on the angles between the pulsar axis of 
rotation, the direction of the beam and the observer pulsars can b e radio-quiet. This source 
of incompleteness can be very relevant: for instance, Ref. 205| estimates the fraction of 
pulsars which are actually beamed in the direction of an observer at Earth, as a function of 
the pulsar period in seconds, as f(P) = 0.09[log(P/s - l] 2 + 0.03. In addition, the Reader 
should bear in mind that the same incompleteness from beaming might be true for other 
wavelengths, including X-ray and gamma-ray emission, although it is matter of debate how 
to estimate the incompleteness in those frequencies, and more specifically how this would 
impact those sources relevant for the present discussion. For sure, the Fermi Large Area 
Telescope (LAT) is giving, and will keep giving in the near future, a crucial contribution 
to the discovery of numerous new bright gamma-ray pulsars (potentially radio-quiet) that 
could very well contribute to the local flux of energetic (see sec. UXj) . 

Fig. [E] shows the positron fraction (left panels) and the total e flux for different models 
for the e pulsar output (top panels), for the injection spectral index a (middle panels) and 
for the diffusion setup (bottom panels). As evident from the discussion in sec. HI1 the ST 
and CCY models give comparable e^ outputs, with the latter boosting the output of more 
nearby pulsars (those contributing to the highest energies). The ZC and the HR model 
predict an almost negligible contribution to both the positron fraction and the overall e flux 
from distant pulsars. In general, in all scenarios, the positron fraction starts to increase after 
100 GeV (unless one assumes a MIN diffusion setup), so it is unlikely that the PAMELA 
positrons dominantly originate from further away than 1 kpc. 

The spiky shape of the curves of fig. fT2l should be interpreted with care: the reader should 
keep in mind that we are assuming a hard spectral index, and no cutoff in the power-law 
injection spectrum at high energies. In addition, we neglect the temporal profile of the 
e ± injection, which again would smooth out the spikes we find. Lastly, the finite energy 
resolution of any e* 1 detector would also make the spikes appear much smoother in the 
actual data. 

The middle panels show that no significant qualitative variation is induced by changing 
the injection spectral index from a = 1.6 to 1.8, 2.0 and 2.2. The harder the spectral 
index, the steeper the spikes predicted in both the positron fraction and in the overall 
differential spectrum. 

A more significant role is played by the diffusion setup. The MIN model drives energetic 
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FIG. 13: Break-up of the contribution from distant (> 1 kpc) pulsars in age (left) and distance 
(right), for a = 1.8 and the MED diffusion setup. 

to lower energies much more effectively than the MAX and MED models, which enhances 
the contribution of distant pulsars to positrons in the PAMELA energy range. The MED 
and MAX setups, on the other hand, give comparable predictions for the contribution of 
distant pulsars. 

In fig. [131 we break up the contribution from distant pulsars to the e spectrum in intervals 
of age (left) and distance (right), assuming a MED diffusion setup, the ST model and a = 1.8. 
As expected, the youngest pulsars contribute to the spike at large energies; pulsars with ages 
between 0.2 and 0.5 Myr contribute dominantly around 1 TeV, while in the range of the 
ATIC bump the bulk of the contribution from distant pulsars stems from those objects with 
an age between 0.5 and 1 Myr, as already observed above. Finally, older pulsars give a 
subdominant contribution at even lower energies. 

The break-up with distance highlights that the ATIC bump is plausibly mostly dominated 
by the nearest objects, with some contribution from pulsars as distant as 2 kpc. Sources 
further than 2 kpc give a negligible contribution, spread out through the entire energy 
spectrum, given the possible age range. 

The bottom line of this section is that the models we used for the computation of the 
pulsars' e ± output are consistent with all relevant objects in the ATNF catalogue (pulsars 
with g > 1 give an even smaller contribution), and distant sources, while unimportant or 
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FIG. 14: The anticipated data on the e spectrum from one year of full science operation with the 
Fermi-LAT telescope, for the two cases, shown in fig. \TT\ of one bright nearby source (Geminga) 
and of multiple, more distant pulsars accounting for both the PAMELA and the ATIC data. 

subdominant for the PAMELA energy range, might indeed be relevant to account for the 
spectral feature observed by ATIC. 



IX. THE ROLE OF FERMI-LAT: SPECTRAL MEASUREMENT AND NEW PUL- 
SARS 



It has been recently pointed out that the Large Area Telescope (LAT) on board the 



Spa ce Telescope can efficiently detect energetic cosmic ray electrons and positrons 



206 



ermi 



207 . 



2081 ] . Being a gamma-ray telescope based on pair-production, the LAT is, intrinsically, 



a high-energy electron and positron detector. The main problem one faces in measuring 
the spectrum with Fermi-LAT is to efficiently separate electrons and positrons from 
all other cosmic-ray species, mainly protons. The LAT team has successfully argued and 
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demonstrated that this is actually possible, with an appropriate choice of analysis cuts. In 
the energy range between 20 GeV and 1 TeV, the effective geometric factor (for electrons), 
after applying the mentioned cuts, ranges from 0.2 to 2 m 2 sr (enormously larger than both 
ATIC's and PAMELA'S, especially when one additionally factors in the data-taking time as 
well), and the energy resolution is 5 to 20%, depending on the energy. Applying the selection 
cuts to the simulated c osmic ray flux , the residual hadron contamination was estimated to 



be « 3% of the e ± flux 



206 



207 



208 



We use this number to evaluate the systematic error 
on the flux measurement in addition to statistical uncertainties. 

Fig. [14] illustrates the capabilities of Fermi-LAT to measure the electron-positron spec- 
trum in one year of full science mode data taking, and assuming a geometric factor of 1 m 2 
sr. The energy bins are logarithmically evenly spaced, but they fall in the ballpark of the 
anticipated detector energy resolution (AE e ±/E e ± ~ 5...20 %). To estimate the error bars 
we linearly sum the systematic uncertainty from hadronic contamination and the statistical 
uncertainty (the former in black, the latter in blue or red, and assumed to be normally 
distributed, and almost invisible at low energies). Systematic errors dominate at low en- 
ergies, while at larger energies the e ± counts get low enough that statistical errors become 
dominant. Both statistical and systematic errors will improve with more data taking time 
and with a better understanding of the detector. The take-away message is, though, that 
Fermi will provide an exquisite measurement of the e ± spectrum. In the figure we use the 
two models of sec. IVHI for illustration. It is clear that in both cases Fermi-LAT will have 
the sensitivity to (indirectly) observe a drop-off in the positron fraction (blue curve), if this 
occurs where predicted in our simple scenario that fits the PAMELA observations, and will 
definitely solidly confirm or refute the spectral feature observed by ATIC. 

Fermi-LAT will play a second decisive role in understanding the origin of energetic e^: 
the discovery of new gamma-ray pulsars. In addition to monitoring portions of the gamma- 
ray sky known to host radio-pulsars, the LAT Collaboration is actively and successfully 
carry ing o ut a program of blind search for gamma-ray pulsars, along the lines outlined in 
Ref. 211]. The first Fermi-LAT science paper was in fact about a newly discovered very 



bright gamma-ray pulsar in the CTA 1 region 



2091 ] . Very accurate measurements of the 



period and of its derivative made it possible to derive a rather acc urate estimate of the age 



of the CTA 1 pulsar, which was found to be T ~ 1.4 x 10 4 yr 



2111 ] . This age determination 



is compatible with the radio and X-ray pulsation of the objects, which also inform us on it 



42 



1 1 1 1 — i — i — i — r~| 

CTA-1, a=1.6, MIN diffusion setup 



r f t- * hi 



100 



> 
a 



X 

m 



10 



Dist- 


dkpc, ST,f=0.1 


Dist- 


■ 1.2 kpc, ST, f=0.3 


Dist- 


=1.4 kpc, ST, f=0.5 


Dist- 


=1.4 kpc, CCY 



_I I I I I I I L 




100 



1000 

Energy [GeV] 



10000 



FIG. 15: Predictions for the contribution to the total flux of the CTA1 pulsar, recently discov- 
ered by Fermi-LAT 209]. We assume the MIN diffusion setu p and a spectral index a = 1.6, as 



inferred from the EGRET data on the same gamma-ray source 



210l | . The thinner blue dashed line 



corresponds to the prediction for a cutoff at 10 TeV in the injection spectrum. 

having a distance of 1.4 ± 0.3 kpc. In addition, the LAT data indicate a rather powerful 
spin-down luminosity of 4.5 x 10 35 erg/s, and a surface magnetic field i? 12 ~ 11. 

Pulsars like CTA 1 are likely to contribute significantly to the e ± flux. To illustrate this 
point, we plot in fig. [15] the predictions for the overall e ± differential spectrum expected 
from the CTA 1 pulsar, for three values of the distance (1, 1.2 and 1.4 kpc) and the ST 
model (with values of f e ± = 0.1, 0.3 and 0.5, respectively), as well as for the CCY model, 
assu ming a distance of 1.4 kpc. Consistently with the EGRET spectrum of the same source 



210| we assume a = 1.6. We employ in this plot a MIN diffusion setup, and we show with 



the thin dashed blue lined the effect of an injection energy cutoff at E e ± ~ 10 TeV. Given 
that CTA 1 is a relatively young object, the contribution peaks at very high energy (at or 
more than 10 TeV). If the source is located between 1 and 1.2 kpc, it plausibly contributes 
to the highest energy bins in the H.E.S.S. data. Fermi LAT has therefore the ability to 
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FIG. 16: Predictions for the anisotropy in the arrival direction of cosmic ray leptons from nearby 
pulsars, defined as (If — lb) /(If + lb), where If is the total number of events from the emisphere in 
the direction of the pulsar, and If, that from the opposite emisphere. The parameters are the same 
as in fig. FH and IT5| blue thin dashed curve for CTA1. For Vela and CTAl we consider a cutoff at 
10 TeV in the injection spectrum. 

discover numerous sources which are bright enough to give a significant contribution to the 
e ± spectrum. 

Besides a very accurate measurement of the high-energy cosmic-ray lepton spectrum and 
the discovery of new local, bright pulsars, Fermi-LAT data will also provide directional 
information on high-emergy electrons and positrons. Although charged cosmic ray leptons 
are isotropize quickly in the ISM magnetic fields, if they are produced close enough to 
Earth, and with large enough energies, a potentially detectable dipolar anisotropy might be 
detected from the direction of nearby pulsars, if indeed the l atte r source the bulk of the local 



high-energy c osmi c ray leptons (on this point see also Ref. 203| and the recent estimates of 



Rcf. 



35| and [l94j |)- 



We show in fig. [TBI the predicted anisotropy (defined as usual as(l/ — i&) / (// + lb), where 
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// is the total number of events from the emisphere centered on the direction of the pulsar, 
and h that of events from the opposite emisphere), as a function of energy, produced by 
the pulsars we considered in sec. II I II and in the present section. Specifically, we adopt here 
the same parameters as those employed in fig. H] for Monogem, Geminga, B0355+54 and for 
Vela (for the latter we assume a cutoff at 10 TeV in the injection spectrum), and, for CTA1, 
the setup we employed for the thin dashed blue line (distance of 1 kpc, ST model with 
f e ± = 0.1, and a cutoff at 10 TeV in the injection spectrum). In the Fermi-LAT range, the 
largest anisotropies lie between one part in a thousand and one in a hundred, for the pulsars 
and the particular setup we employ here. Notice that the distance of B0355+54 suppresses 
the resulting anisotropy as opposed to a nearby object like Monogem, when comparing fig. [TH] 
with the fluxes shown in fig. HI 

No dedicated studies of the performance of Fermi-LAT in discerning an anisotropy in the 
arrival direction of high-energy cosmic-ray leptons exist yet, so it is hard to estimate whether 
anisotropies at this level might in fact be detected. Naively, one can conservatively estimate 
the number of events from the diffuse cosmic ray lepton background alone as a function of 
energy, and compare, in a given emisphere, a 2 or 5cr poissonian fluctuation in that number 
to the anisotropy to be measured. We show curves, obtained in this way, corresponding to 
1 yr of data, and to 2 and 5a, as well as two 2a curves for 5 and 10 yr. 

The best chances Fermi-LAT appears to have to detect one of the bright nearby pulsars is 
(i) for a mature and nearby object like Geminga, in the low energy part of the spectrum; (ii) 
for a younger, but farther bright pulsar like Monogem, in the several hundred GeV range, 
or (iii) for very young pulsars, like Vela or CTAl, in the mufti- TeV range, if the systematics 
can be controlled enough to reconstruct with accuracy electron events at those energies. 
This latter case can be in principle nicely complemented with data from ACT's as well. In 
any case, it appears that an anisotropy study of the Fermi-LAT is potentially feasible and 
even promising. 

We point out, however, that while the detection of an excess from the direction of a 
nearby pulsar would be a very suggestive experimental result, the absence of an anisotropy 
feature at the level predicted here would not be a conclusive evidence against the pulsar 
scenario (for instance due to local, solar system effects, or to the large scale structure of 
ISM magnetic fields, or because of the proper motion of the pulsars themselves, and thus 
uncertainties in the actual site of e ± injection). 
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To conclude, we wish to stress here that the Fermi telescope will play an essential role 
in forming a self- consistent picture of local high-energy cosmic ray electrons and positrons, 
specifically in connection with, but not limited to, the pulsar scenario. 

X. SUMMARY AND CONCLUSIONS 

We argued in this study that there is indeed no "need to add entities beyond necessity" 
to explain the PAMELA and the ATIC data on cosmic-ray electron and positron fluxes: 
existing and well-known pulsars easily and naturally (from an energy budget standpoint) 
account for both experimental results. 

We summarize below the main results of the present theoretical study. 

• Local pulsars are well-known sites of electron-positron pair production. Although a 
firm theoretical understanding of the pulsars' e ± energy output is lacking, we surveyed 
here a few possibilities, which are all approximately in line with one another and with 
estimates of the output energy in electrons from SNR's. Remarkably, by picking a few 
nearby known pulsars and employing the mentioned energy output models, we easily 
reproduce the spectral features and intensities reported by both PAMELA and ATIC, 
with nominal values for the e ± output efficiency; 

• The quality of the recent data on the positron fraction and on the differential e ± flux 
is good enough to allow the exploration of the regions of the pulsar parameter space 
favored by the two experimental observations; uncertainties in the determination of 
pulsars' ages and distances - particularly, as far as the latter are concerned, stemming 
from the large observed pulsar velocities - can however significantly compromise our 
ability to pinpoint single major contributors to the local cosmic-ray lepton flux 

• Although statistically and energetically plausible for existing objects, we find that the 
PAMELA and ATIC data are not well fitted by one single bright nearby source; 

• For a given e^ injection spectrum and diffusion setup, the PAMELA data point towards 
a corridor of correlated values of pulsar age and distance. Several known pulsars fall 
within this corridor, and could coherently add up to contribute to the positron excess 
over the standard secondary population observed by PAMELA; 
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• The ATIC feature at around 600 GeV can also be explained by one or more existing 
pulsars, although this generically requires a rather high electron-positron energy out- 
put, and, unlike PAMELA, it does not seem to be a natural or necessary post-diction 
of the pulsar contribution to energetic e^: while it is plausible to expect pulsars to 
produce a positron fraction like that measured by PAMELA, it is not as obvious, in the 
same context, to expect a spectral feature like that reported by ATIC. Nevertheless, 
relatively young (0.4-0.6 Myr) pulsars 1-2 kpc away can easily account for a spectral 
feature such as the one observed by ATIC. 

• The Fermi Space Telescope will play an extraordinarily important role in the under- 
standing of the local flux of energetic electrons and positrons, by (1) providing us, in 
the very near future, with an exquisite measurement of the e ± flux, that could unveil 
several interesting features, and/or rule out recent experimental claims, and (2) by 
discovering gamma-ray pulsars that can add-up to give a self-consistent picture of the 
origin of the local energetic electrons and positrons, and (3) searching for anisotropics 
in the arrival direction of high-energy electrons and positrons. 

A convincing picture of the nature of energetic cosmic ray electrons and positrons evidently 
calls for a two-pronged effort: on the one hand, theorists working on cosmic rays and on 
compact objects should refine the theoretical understanding of production in pulsars 
and possibly in other astrophysical galactic objects; on the other hand, further and deeper 
observations of candidate electron-positron sources at various wavelengths (including radio, 
X-ray, gamma-ray and very high energy gamma-ray frequencies) will be needed to comple- 
ment and to drive the theoretical effort towards the construction of a thoroughly satisfactory 
picture. 
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